A step towards testing general relativity using weak gravitational lensing and redshift 

surveys 



o 
o 

(N 

< 
(N 



I 

O 



> 

O 
O 

o 

oo 

o 



X 



Yong-Seon Song^ and Olivier Dore^B 

^Institute of Cosmological & Gravitation, University of Portsmouth, Portsmouth, POl 3FX, UK ana 
^ Canadian Institute for Theoretical Astrophysics, 60 St. George St, 
University of Toronto, Toronto, ON, Canada M5S3H8 
(Dated: August 24, 2009) 

Using the linear theory of perturbations in General Relativity, we express a set of consistency 
relations that can be observationally tested with current and future large scale structure surveys. 
We then outline a stringent model-independent program to test gravity on cosmological scales. 
We illustrate the feasibility of such a program by jointly using several observables like peculiar 
velocities, galaxy clustering and weak gravitational lensing. After addressing possible observational 
or astrophysical caveats like galaxy bias and redshift uncertainties, we forecast in particular how well 
one can predict the lensing signal from a cosmic shear survey using an over-lapping galaxy survey. 
We finally discuss the specific physics probed this way and illustrate how f{R) gravity models would 
fail such a test. 

PACS numbers: 



I. INTRODUCTION 



To understand the origin of the accelerating expan- 
sion of the universe is one of the sahent question in con- 
temporary cosmology [E S H, 0, H, H, 0|- com- 
monly attributed to the existence of some extra unknown 
physics loosely labelled Dark Energy (DE) Q. Alter- 
natively, it might cast some doubts on the theoretical 
foundation of current cosmology, that is the theory of 
General Relativity (GR). From a theorist point of view, 
as most efforts to construct self-consistent DE models 
within GR seem unsuccessful, there is no a priori rea- 
son to rule out a modification to gravity on cosmologi- 
cal scales as an explanation [9]. From an observer point 
of view, whereas GR passes direct tests probing the So- 
lar System scales (lO^^m) down to the laboratory scales 
(10~'^m) jlOj], testing gravity on cosmological scales is 
challenging and overly model dependant so far. In this 
work, we propose to remedy these issues by defining a 
set of simple model independent tests of GR on cosmo- 
logical scales, although more complicated models of GR 
(clumping dark energy or interacting dark energy mod- 
els) are not considered in this paper (the complete set 
of test including both exotic dark energy models is be- 
ing prepared as the second step of testing gravity) . The 
tests we devise rely on simple and testable consistency 
relations based on linear GR and make use of the appro- 
priate combinations of various observables of large scale 
structures. Were any of these relations proved to be vi- 
olated, it would be a clear sign of a break-down of GR 
and it would highly motivate theoretical work on alterna- 
tives to GR. Reversely, before any such conclusion can be 
reached, it requires an absolute demonstration that ob- 
servational systematic effects are well under controlled. 



A firm believer in GR might look at these states as sim- 
ple test of systematic effects. Note that throughout this 
work, we will assume the validity of standard cosmology, 
that is to say the validity of the Copernican principle [ll| . 
If this was not the case, then the DE problem could be 
solve without new physics and more specific tests would 
have to be designed [l^, [l^, ■ 

Several alternative approaches to test GR have been 
proposed so far, some parametric, some non-parametric. 
The parametric approaches can be classified into sev- 
eral categories according to their goal: i) to separate 
the geometrical and growth signatures in the w <k Wa 
plane [l^, [l^ [13], ii) to define the probed 'w{z) range 
in a model independant manner [isl . [l9l [20| , iii) to look 
for an anomalous linear growth rate using the 7 parame- 
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iv) to parametrize the 
M 132, 33, 34, 35, 3^, 



ter [191, 1211, 123,1231, 12; 

metric perturbations 

I37I , [ssl [39l . l40j . Alternatively, non-parametric approaches 
were developed by [Sllli^ |43|. l4i| . 



In this paper, we advocate the use of a non-parametric 
approach - based on a set of GR based consistency rela- 
tions - to test gravity on cosmological scales. We will 
use explicit alternative to GR for illustrative purpose 
only. To do so, we jointly use various observables, i.e. ve- 
locities, galaxy clustering and weak-gravitational lens- 
ing [4l|. More precisely, by using a combination of the 
two first, we will be able to compare it to the later. In 
this work, extending the work of [42i | we focus closely on 
potential observational biases. We propose in particular 
various ways to deal with the galaxy bias, redshift space 
distortion and spectroscopic and photometric redshifts. 



We first begin by introducing in SecITTla general frame- 
work to test GR on cosmological scales. We then detail 

how to build faithful probes of matter perturbations in 

Sec. Uni before detailing one test of GR in SeclIVI We 

* Electronic address: |Yong-seon.song@port.ac.uk,olivier@cita.utoronto.ce[ 3nclude and discuSS our findings in Scc. IVl 
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II. GRAVITATIONAL CONSISTENCY TESTS 
ON COSMOLOGICAL SCALES 

In a metric theory of gravity, large scale structures cor- 
respond to metric and matter-energy momentum pertur- 
bations. In the Newtonian gauge, metric perturbations 
are described by 



ds^ = -(1 + 2^')dt^ + (1 + 2<!>)a'^Sijdx^dx^ 



(1) 



where $ and 4* denote curvature perturbations and the 
Newtonian force respectively. In this work (except in 
section Sec. IIVBI) we will aim at constructing GR con- 
sistency relations and thus we will work within GR. We 
will restrict ourselves to a linear theory of perturbations 
on sub-horizon scales to get simpler and easier to test re- 
lations. We will also consider matter density fluctuations 
to be dominating the growth of perturbations and thus 
neglect DE clustering. 

The first simple relation that can be tested within GR 
is based on the continuity equation: 



dt 



(2) 



where Om is given by 0„j = V • Vm- Both sides of this 
equation can be probed observationally independently 
and we shall call the observational test of this relation 
the energy-momentum consistency test. As will be seen 
below, whereas peculiar velocities can be traced directly, 
probing matter density fluctuations requires to deal with 
several observational artifacts. In this paper, we will as- 
sume this relation to be satisfied and we will use it for 
our observational determination of the galaxy bias. Note 
however that it could be violated due to dark sector in- 
teractions. 

The second consistency relation stems from the lack of 
anisotropic stress and relates metric perturbations. Since 
in this paper, we neglect DE clustering, only the matter 
component clusters and the no anisotropic stress approx- 
imation is valid, 



$ + * = 0, 



(3) 



which reduces the degrees of freedom of metric perturba- 
tions. While the Newtonian force 5* can be reconstructed 
from the evolution of peculiar velocities, the curvature 
perturbation $ is given by matter fluctuations. Thus 
both observables could be compared to determine the 
presence of non-trivial anisotropic stress which is pre- 
dicted in most modified gravity models and dark energy 
clumping model. We call this test the metric consistency 
test. 

The other tests relate metric perturbations to matter- 
energy fluctuations, dynamically or non-dynamically. 
Newtonian force sources the dynamics of matter fluctu- 
ations. If the time evolution of peculiar velocity can be 
reconstructed, then so can '5 through the Euler equation. 



dOm 
~dF 



-HOn 



(4) 
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FIG. 1: Web of cosmological tests of GR (see an analogous 
plot in Fig. 1 of ;3^J). 



which we label as the dynamical constraint test. Because 
the degrees of freedom available in metric perturbations 
are reduced by the no-anisotropy condition in GR, we 
shall not use this constraint in this paper. 

Finally, the relation between curvature perturbations 
and matter fluctuations yield the non-dynamical con- 
straint test, also known as the Poisson equation. It is 
a key relation which we will use and test in this work. It 
writes as. 



(5) 



Using previously defined relations, in principle both sides 
of Eq. Oare calibrated and can be readily used to test it. 

This set of relation constitute a web of possible cosmo- 
logical as illustrated in Fig. [1] 

As we hinted at before, to test any of this relations, 
we need various observables. To trace matter pertur- 
bations, Smj we will use biased tracers like galaxies or 
clusters, whose relevant observational caveats will be ad- 
dressed below. Velocity surveys, direct (kinetic Sunyaev- 
Zel'dovich or peculiar velocities measurements) or indi- 
rect (galaxy redshift survey) allow us in principle to probe 
Qm- Finally, metric perturbations are directly probed 
by weak gravitational lensing which distorts the shape of 
source galaxies along the line of sight. Since the geodesies 
are determined by the gradient of combination of both $ 
and as V(<I' — ^'), weak lensing will probe the inte- 
grated effect of metric combination which is 2V$. We 
will introduce our notations for those observables below 
before discussing in more details the specifics of the test 
we are interested in. 
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III. TRACING MATTER PERTURBATIONS 
WITH GALAXIES 

Whereas there is no direct probe of matter density fluc- 
tuations (even though pecuhar velocity measurements 
are usually considered to be a direct tracer of density 
fluctuations, this requires the energy momentum con- 
servation law written in Eq. [2] which is not granted in 
our context), galaxies can still be used for this purpose 
provided several observational artifacts are properly ac- 
counted for. 

Since we are working on large linear scales, we will ig- 
nore non-linear effects which usually are another limiter. 
We will nevertheless discuss in the next sections redshift 
distortions, redshift measurement uncertainties and bias. 
The information extracted from galaxy density fluctua- 
tions will be limited by how well we can handle those 
effects. Because the tests we are presenting require the 
use of projected quantities, we will focus in particular 
on how those artifacts affect the projected galaxy dis- 
tribution within a redshift bin in its comparison to the 
projected matter distribution. 

In this paper we will restrict ourselves to compare 
the measured and predicted power spectra, even though 
a map-based approach would be feasible. As we will 
work mostly on large scales where the field distribution 
is Gaussian, considering this approach is lose-less. It is 
certainly not true on smaller scales. 



A. Cross-Power spectra conventions 

In this work, we are interested in cross-correlating vari- 
ous tracers of the curvature perturbations at various red- 
shifts to test gravity. We will consider three observables 
based on gravitational lensing, galaxy counts and velocity 
measurements from rcdshift-space distortion. In particu- 
lar, we will sometimes combine some projection of those 
observables and we will thus define the associated angu- 
lar power spectra below. 

We will consider that the Limber approximation is 
valid for all those observables in the angular range of 
interest. The projected angular power spectra Cf-^ of 
any pair of perturbations X and X' is given by 



tion n(z), defined as 

ng{z) oc z2e-(-/i-5)^ 



(7) 



XX' 



25 (£-f 1/2)3 
X J dDDW''{D)W'''{D)Aiciao,k^ -^){6) 

where A^(;(ao,/c) is the rms amplitude of curvature 
fluctuations on comoving hyper-surfaces at the starting 
epoch oo during matter domination and D is the angular 
diameter distance. The window function (D) is de- 
termined by the property of the quantity X as is shown 
explicitly below. 

We consider that galaxies follow a continuous distribu- 



where we assume the space-based surveys [43]. We con- 
sider that for each observed galaxy, we can measure either 
its spectroscopic redshift or its photometric redshift. We 
neglect the errors on the former and we will discuss the 
latter further in Sec. IIIID2I 

The galaxy over-density, 6g, are measured on the red- 
shift slice labeled by i at the comoving distance Di from 
the observer and is assumed to be a linearly biased tracer 
of the matter over-density Sm, with bias b, so that the 
window function for Si^g is given by 

ly^nA) = 3G.-^K,fc)^n,(z.)6(z,) ^24j^ (8) 

where the growth function Gs^ is given by G^^ {ai , k) = 
a<^{a„k)/<^{ao,k). 

The deflection angle d due to gravitational lensing is 
defined by the gradient field of the lensing potential, d — 
V0, where 



(z) = -2 dD 



D(z') - D(z) 



The window function for is thus 

W'^{D{z)) = ^2G^{a,k) dz'n,{z')^^^^^^^J^^O) 

where the growth function G'$(a, k) is given by 

G<i.(a, k) = $(a, A:)/$(ao, k) . (11) 
The angular power spectra of the deflection angle is given 

by c;'''' = /(? + i)c;**. 

Finally, peculiar velocity in our context will be mea- 
sured following the method proposed by Song and Per- 
cival |4^]. On large angular scales, the evolution of LSS 
measured by peculiar velocity is an independent tracer of 
the history of LSS. If we consider the energy momentum 
conservation law written in Eq. [5] to be valid, then the 
window function for 9i „ defined in the i*'' shell is given 

by 

V^-(A)^|c..,(aaO§n,(«.)|±||. (12) 

where Gg^ = aGs^ (^i i k) is the growth factor for 6*^ and 
the velocity bias is ignored, 9i^g = Oi^m- 



B. Redshift distortion and angular power spectra 

In this section, we consider a spectroscopic survey 
whose redshift measurement errors are neglected and dis- 
cuss the effects of redshift space distortion on the angular 
power spectra. The galaxy density fluctuations measured 
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in redshift space are distorted by peculiar velocities on 
all scales. The observed power spectra can thus be a 
mixture of density fluctuations and peculiar velocities. 
To derive the projected galaxy angular power spectra on 
large scales, i.e. linear, we first write the galaxy density 
in redshift space as [43] 



5g{k, D) = a,)Wg{D){l + //i') 



(13) 



where curvature perturbations are separated into a scale 
and a time dependant part. 



$(fc, a) = ci>„(fc)l(^ ^ ^o{k)aGsJa) , 



(14) 



and where ap is initial epoch at matter domination and 
Gs^ is growth function of matter fluctuations and the 
window function Wg is defined in Eq. [5) The fractional 
weight function between density fluctuations and peculiar 



velocities / can be written as 

d In Ga,^ (a) 



din a 



(15) 



The projected 2D galaxy density from redshift space 
is then written as 

5g{h) = J dDWg{D){l + ffi^)6,^{n,D) (16) 



We can then write 



''Im 



dn[n)5g{h)Y;^{n) 



(17) 



hence 



(2^+l)Cf = ^(aLa: 



a* \ 

tml 



(18) 



dMf7'dMZ?'^^($(k)<I>*(k'))T^,(l + + /M")e^^'=''^-''''^V;„(n)r,™(n') 



E 



(2^ 



dmDWg{l + ft?)e^'^^''Y;^{n 



J dn'dD'W;^il + /M'2)e^'='''^V,„(n') 



(19) 



We can then calculate the bracket as 



the 80 modes and we have 



[...] = JdndDWg{l + fn')e^'''^^Y;^{h) 

= v/7r(2« + l)5mO J (iWDl¥<,(l + V)e''^^P,(/z) 



2v/^(2Z + l)S,nO / dDWgiD) 



Ji{kD) - f ( ji-2ikD) - —ji_,{kD) 



We thus obtain 



where 

m 



= 47r / ^p^{k)[mf 



dDWgiD) 



.(20) 



(21) 



(22) 



JiikD) - / ( ji-2{kD) - —ji^kD) 



At high £ modes where the Limber approximation is safe, 
we have kD ~ £ which suppresses the contribution from 



/z(fc)~ / dDWgiD) ji{kD) , 



(23) 



where the redshift distortion effect on the projected 
power spectra is suppressed. 

Those two equations illustrated in Fig.[2]show that pro- 
jecting the angular power spectrum minimizes the red- 
shift space distortion through cancellation effects. This 
effect is interesting here since it makes the projected an- 
gular density a more direct tracer of the projected matter 
density. This will not be valid anymore when non-linear 
effects are non-negligible, that is for scales larger than 
k - 0.3 Mpc"^ 



C. Photometric uncertainties and angular power 
spectra 

We now consider a photometric survey and discuss the 
effects of redshift errors on the projected galaxy angular 
power spectrum. We assume for the latter Gaussian error 
with rms 



(t(z) = aphz(l + z) 



(24) 
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FIG. 2: Projected galaxy over-density angular power spec- 
trum as a function of multipole £ as defined in Eq. [21] and 
Eq. [23] (solid line) and its limber approximation as defined in 
Eq. [23] (dashed line). The various curves correspond to vari- 
ous bin width. Obviously, the wider the redshift bin the better 
is the Limber approximation but it is not a very accurate one 
except at smaller scales. 



A subset i of galaxies whose photometric redshift is such 
that <2< Zi, follows the distribution 



variance are widely different. Since the curvature pertur- 
bation power spectrum is weighted by the square of ng{z) 
spectrum when Eq. [6] is applied to the galaxy overden- 
sity, 5g , this introduces the substantial bias we see in this 
plot. If the true underlying distribution was known, i.e. 
if wc could deconvolve the observed distribution with the 
photometric error distribution, then we could device eas- 
ily weights that do not lead to such an error. In practice 
though, it is unlikely that this deconvolution will be fea- 
sible and we thus resort to another way to correct for this 
effect by introducing another bias factor that we called a 
photo-0 bias, hz- 

Since Ug is a smooth function of z and since the mat- 
ter power spectrum is mostly featureless, we expect bz 
to be weakly dependent on scales. In the right panel of 
Fig. [31 we plot the relative difference between the true 
galaxy angular power spectrum, Cf^'^^^ , and Cf^ with 
Uz — 0.03. In the low redshift bins, there arc non-trivial 
scale dependent in the angular range coming from bary- 
onic features (lOOMpc at z = 0.3 corresponds roughly to 
t = 60). In the high redshift bins, the scale dependence is 
ignorable. Fortunately, the contribution from those low 
redshift bins to the reconstructed lensing power spectra 
is not significant. It means that we can treat in a first 
approximation the biasing due to the limited photo-z er- 
ror as scale independent. Within this hypothesis, Cf^ is 
a linearly bias tracer of the matter angular power spec- 
trum, with a total bias = ^g^z; where bg denotes the 
bias due to galaxy and bz denote the bias due to photo-z 
uncertainty. The measure of the 3-D galaxy power spec- 
trum allows us to measure bg separately but we can also 
measure directly the total bias in a redshift bin, i, as will 
be detailed in Sec. IIIID2I 



n.i{z) 



A,. 



ng(z) 



erfc 



Zj-l - z 

V2a{z) 



— erfc 



/ Zj- z 

\V2<j{z) 



where erfc is the complementary error function and Ai 
is determined by a normalisation constraint. Through- 
out this work, we will choose Ai so that J dzni{z) = 1. 
Doing so, we discard the total number of galaxies as an 
observable. 

As an illustration, we consider three different levels of 
photometric errors, respectively cr^ = 3 x 10"'', 3 x 10"'^, 
3 X 10~^. While the first noise level seems idealistic, 
the second one seems achievable in a near future (e.g. 
[49| ) and the third one corresponds to what is currently 
achieved with SDSS [1^. The corresponding window 
function ng{z) defined in Eq. [^IJis displayed in the left- 
top panel of Fig. [3] The left-bottom panel of Fig. [3] shows 
the simple projected angular power spectra of the galaxy 
overdensity according to Eq. [6] As we can see, the pres- 
ence of photometric errors introduces an important error 
that has to be taken into account. 

This error can be simply understood the following way. 
Looking at the left panel of Fig. [3] one sees that whereas 
the mean of each distribution and their integral - the to- 
tal number of galaxies in each bin - are identical, their 



D. Galaxy bias and associated uncertainties 

We now present two alternative ways to measure the 
linear galaxy bias, assuming either a precise spectroscopic 
survey or a weak-lensing based photometric survey. 



1. Using a spectroscopic survey 

Using precise spectroscopic redshift measurements, we 
are able to separate the peculiar velocity power spec- 
tra from the redshift-space power spectrum P°'^^(k) of 
a galaxy redshift survey [40|. The latter is commonly 
modeled as 



Pf^(k) 



Pggik) + 2M2Pge,(k) + /i'e.e.lk) 
X F{eal{z)fi^). (25) 



The separation of is Pgg(k) and Pee(k) made possible 
using the angular dependence of Pg^^ (k) and where O = 
9/aH 51]. 

We quantify how well this separation can be performed 
using a Fisher matrix formalism. The Fisher matrix anal- 
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FIG. 3; Left panel, top: Ug distribution plotted for three various a^, i.e. 3 x 10'^, 3 x 10"^, 3 x 10~^ and the underlying galaxy 
distribution. Bottom : Corresponding galaxy overdensity angular power spectrum, C®". Photometric redshift error entails 
an important bias. Right panel: Relative diflterence between Cf^ power spectra with or without photo-z errors considering 
(Jz = 0.03. Solid, long dash, short dash and dotted curves correspond respectively to Cf^ at 2=0.3, 0.5, 1.5 and 2.5. 



ysis relevant to this separation for a given k and z bin is 
given by [Slj 



k^dk 



2(2^)2 
g v^,f^,Zj) dlnP° 



d In P°'^"(fc, fi, Zj) d In P°^''{k, zj) 



dpa 



dp0 



dn 14ff(fc,/i, Zj) 



(26) 



where a and /3 run from 1 to 2 and denote respectively 
Pgg and Pee and Wp^Q is a down-weight function filter- 
ing out the modes contaminated by non-linear redshift 
space distortions. Note that PgO^ is considered to be 1 
in this separation which is valid on the linear scales of 
interest to us. The effective volume V^g in each redshift 
bin j is 



Vesih, ^J,,Zj) 



.^P°^^{k,,fl,Z,) 



-\ 2 



niP°^^{k„fi,Zj 



1 



survey 



{z^ 



(27) 

where is the shot noise term coming from the finite 
galaxy density supposed constant here, and Vsurvcy (zj) 
is the survey volume in a given redshift bin. For the 
large scales of interest to us, the cosmic variance term 
dominates over the shot noise and Vesiki, zj) is nearly 
identical to T^urvey (-Zj)- For our estimation, we will con- 
sider a full sky survey with a constant galaxy density of 
n = 5 X 10~^ h^Mpc~^ and a constant bias equals to 



1. As illustrated in left panels of Fig. [4l we are able to 
separate properly Pgg and Pee for wide k and z bins. 

Following this measurement of Pgg and Pee, we can 
constrain simultaneously the cosmological parameters 
and the galaxy biases. If the galaxy bias is scale- 
independent and depends only on redshift, then the 
real space galaxy power spectrum Pgg can also be writ- 
ten in terms of the fundamental cosmological parame- 
ters plus a vector of bias parameters. If we consider 
for example a survey up to z = 3.2 with 8 redshift 
bins of width Az — 0.4, the standard DE cosmologi- 
cal parameter set - as the one used in the Dark En- 
ergy Task Force report [s^l - is extended to 16 elements, 
q = {w,Wa,Wm,Wb, As,ns, 0S,b-i=i-8). For this 

extended cosmological space, the Fisher matrix simply 
writes 



Frr 



ij al3 



dqm 



(28) 



Since Pee is independent of bias, the degeneracy between 
cosmological parameters and the bias is broken [ssj , and 
we can simultaneously measure in each redshift bin the 
bias and e.g. the dark energy parameters. The resulting 
bias uncertainties are given for each redshift bin in Ta- 
ble 1111 D II and are typically at the percent level. Note 
that the galaxy bias is measurement is not detrimental 
to the dark energy parameters as illustrated in the right 
panels of Fig. U) 
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FIG. 4: Left panels: The top panel show the fractional errors for the reconstructed Pq q using the Fisher matrix formalism written in 
Eg. 1261 The bottom plot shows the corresponding fractional errors for Pgg- Right panels: 68% CL contour plots in the w — Wa plane. 
The top panel shows the constraints obtained using Pq^q^ (no bias marginalization) only and Pq q +^99 (with bias marginalization) . 
The bottom panel shows the constraints using Pg, q ~^Pgg ^f^d various value of the parameters ath that quantifies the accuracy of the 
modeling of the Finger of God effect (see Eq. 15 of White et al. 51] for details). 



TABLE L The fractional error of bias in some selected redshift 
bins using Eq |28l 
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2. 


Using 


a photometric survey 







Some promising SNR numbers are given in Tab. HIl where 
we used all the £ up to £ = 500, assumed = 0.03 
and considered a redshift binning for the WL signal from 
z=0 to 3.2 with Az — 0.4. Note that for this evaluation, 
unlike in the previous section we did not vary the other 
cosmological parameters but here again, we expect the 
DE parameters to be non-degenerate with the bias when 
we include the the projected galaxy and weak gravita- 
tional lensing cross-correlation. 



As introduced earlier, redshift uncertainties as big as 
the ones resulting from photometric redshift measure- 
ments introduce an extra, almost linear bias, 6z- Since 
this bias does not affect the cross-correlation between 
galaxy and peculiar velocity, we can only determine the 
"total bias" by cross-correlating the weak gravitational 
lensing signal and the projected galaxy density. The 
signal-to-noise ratio for this correlation is simply given 
by in 



8 /mB, 

EE 



/sky(2Z+l)(C;^'=)2 



and the fractional error on the galaxy bias is 

Abt ^ J_ 
bt ~ S/N ■ 



(29) 



(30) 



TABLE IL Signal to noise ratio estimate for the total bias as 
defined in Eq. [29] in selected bins. 
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IV. CONSISTENCY TESTS 

Now that we have presented how to obtain accurate 
estimates of the projected matter angular power spec- 
trum using galaxy surveys, we proceed to the core of our 
study, that is the details of our cosmological consistency 
tests. We will propose two tests. Either we predict the 
lensing convergence power spectra using a galaxy survey 
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FIG. 5: Relative uncertainties for the reconstructed Cf"^. The dash 
line corresponds to the statistical uncertainties when measuring the 
bias using a large spectroscopic survey (as in Sec. 1111 D 11 1 and the 
dotted line corresponds to the statistical uncertainties when using 
a photometric survey (as in Sec. llll D 2|l . Other systematic bias are 
illustrated in Fig. |6] 



and compare it to the measured lensing power spectra, 
or we predict the cross-correlation between matter and 
galaxy. While the first constitutes an observational im- 
plementation of both the metric test written in Eq. [3]and 
the non-dynamical constrain test, the second one is a di- 
rect implementation of the non- dynamical constrain test 
written in Eq. O 



A. Predicting the lensing power spectra 

The density perturbations are measured on the redshift 
shell labeled by i at the comoving distance Di from the 
observer. In the approximation of a quasi-static evolu- 
tion of perturbations, i.e. considering the perturbations 
constant within a redshift shell, the projected angular 
power spectrum can be written as 



(ji gg _ 



27r^ 



G + 1/2)3 



(31) 

Similarly, the weak lensing power spectra can be dis- 
cretized as 

C!'' = ^ E AAA^%^A**(a., k) .(32) 

If we first assume that there is no dark energy pertur- 
bations, then A$$(ai,A:) can be written in terms of the 



angular power spectrum of galaxies within this shell as 
9 f dz " 



87r2(/ + l/2) AA \dD 



-Uib, 



02 F74 
T — 



(33) 



Substituting this into Eq. ([32]), we are lead to define the 
reconstructed lensing power spectra 
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Note that this estimator is the simplest we can device 
and that we considered the noise to be negligible. We 
also ignore correlations within redshift bins, which is true 
if they are wide enough. If these hypothesis are not full- 
filed, it is straightforward to generalize our estimator to 
handle those effects in an optimal manner. In Fig. [5] we 
show the statistical errors for the reconstructed power 
spectra using our nominal survey. Obviously, the recon- 
struction performs very well. In the right panel of Fig. [BJ 
we plotted several reconstructed power spectra, before 
and after photo-z bias reconstruction for Az = 0.4 bins 
and photo-z errors defined by = 0.03. Obviously, the 
reconstructed estimator agrees well with the input ones 
once corrected from the photo-z bias. As expected fol- 
lowing the results of Fig. [3] though, this bias is harder to 
correct at low z. 

Once this estimator is defined, we can calculate the 
variance of ACf '^'^ as 



which gives a fractional error 

I \ -pi ( 0^hi.\ 

dcT ~ 



1/2 



(35) 
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(36) 



We show the predicted lensing signal accuracy in Fig. [6] 
using the same survey parameters as for Fig. [J] The 
resulting statistical uncertainties in the predicted angu- 
lar power spectrum (right panel) are much smaller than 
the potential reconstruction biases from non-linear effect, 
limited photometry measurement accuracy and discrete- 
ness effects. Percent accuracy is still possible with com- 
ing surveys and we will discuss further in Sec. |V] the sub- 
sequent cosmological interpretation. 

The number of galaxy redshift bins used to approxi- 
mate the continuous kernel of lensing potential along the 
line of sight in Eq. [34] is limited. The choice of the opti- 
mal bin width is motivated by several issues. We chose 
the thickness of the bins to be larger than any correla- 
tion length between two subsequent bins. Besides, since 




FIG. 6: Left panel: Cf'' for a source distribution at 5=0.2, 1.0 and 3.0 (from bottom to top) with a redshift width Az = 0.4. (top 
panel) The solid curves represent the original WL power spectra, and the dash curves represent the reconstructed one when the projected 
galaxy angular power spectra have been cut below some non-linear threshold chosen here to be fc ^ 0. 3 Mpc~^ .{bottom panel) The solid 
curves represent the original WL power spectra, the dotted curves represent the reconstructed ones with photo- z error of cr^ =0.03 and 
the dash curves represent the predicted ones corrected for the photo-z bias. The agreement is obviously good except at low z were the 
photo-z bias is harder to correct for. Right panel: Expected relative uncertainties in the lensing predicted angular power spectra for the 
source distribution at two various redshifts (1.0 and 3.0 respectively from top to bottom panels and again Az = 0.4). The thin curves 
represent the reconstructed lensing signal using all the information available from the projected galaxy density. The visible noise at low 
£ originates in the finite number of galaxy redshift bins available. The thick curves is obtained when removing all the galaxy information 
above k ~ 0.3Mpc~^ to stay in the linear regime. The black curves correspond to the reconstruction for a ACDM model. The blue curves 
to the reconstruction for f{R) gravity with Bq = 1 (see Sec. |V] for more details). Whereas the reconstruction performs really well in the 



the projected galaxy angular power spectra are defined in 
redshift space, a larger width smooth the redshift distor- 
tion effect as shown in Fig. O We find a width of roughly 
Az = 0.1 to be satisfying. The thin curves in the right 
panel of Fig.[6]shows that whereas the accuracy of Cf'^ at 
z = 0.2 is limited by discreteness effect (we can use only 
4 z bins), it is nearly ignorable at high redshift lensing 
bins. 

Another important bias factor involves non-linear ef- 
fects. Whereas the Poisson equation we used to predict 
the lensing signal is valid on all scales, to reconstruct 
a reliable projected density template from the galaxies 
measured in redshift space might be challenging in the 
non-linear regime. To illustrate how important this effect 
is, we filter out the galaxy template for k < 0.3h/Mpc 
and then project the galaxy density field. In the right 
panel of Fig. [SI the stepwise curves represent the result- 
ing predicted lensing signal. Because non-linearities are 
stronger at low redshift, this bias is more important when 
predicting the lensing signal at lower redshift. 

The uncertainty due to photometry measurement can 
be another significant bias as shown in the left-bottom 
panel of Fig. [SI However, we presented in Sec. lIII lT2l how 



to deal with such an effect using spectroscopic surveys 
or photometric surveys. This bias should thus not the 
practical applications of our test. 



B. Adding cross-correlation with velocity 



In addition to the comparison between the predicted 
and the measured lensing power spectra at various red- 
shifts, we can device another consistency test using cross- 
correlations between galaxy and weak lensing. This test 
will highlight in particular any deviation from the Pois- 
son equation that we parametrize with the a parameter 
defined by 



fc^<I> = 4:TrGNCta^ PmSn 



(37) 
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FIG. 7: Black solid curve represent ai for ACDM and blue solid 
curve represent for f(R) theories with Bq = 1. Dash curves 
are errors estimated when using WL-galaxy correlations (as in 
Sec. lIIID2l l and dotted curves are error estimated using the galaxy- 
peculiar velocity correlations (as in Sec. IIII D 



a = 1 corresponds to GR. An estimator for a can be 
simply derived using C; ' and C7 that we rewrite as 
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The estimator a^ is given by, 



dz 
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As such, the fractional error for is equivalent to the 
fractional error for bi given in Tab. IIII D 11 
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A6^ 
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(40) 



and also we can also relate it directly to the compare it 
with the fractional error of '''' using Eg 1551 

In Fig. [71 the estimated a from galaxy maps and cross 
correlation with lensing potential is presented. Solid 
curve is the estimated a which departs from unity around 
edges of lensing kernel sourced by galaxies at z = 3.0. 



But the estimation of a is very close to what is expected 
from GR model, unity. Due to diverse bias effect, it is 
not identical to unity, but it is just off by a percent- 
age level. Still the difference due to systematic effect is 
much smaller than the theoretical deviation due to vio- 
lated GR assumptions, e.g f{R) (R:Ricci scalar) gravity 
model with Bq ^ 1 (long dash curve in Fig. [7]), where Bq 
is defined as, 



d^f/dR^ dR fdlnH 



1 + df / dR din a \ d In a 



(41) 



Dash and dotted curves are statistical errors from dif- 
ferent galaxy bias estimation from redshift survey alone 
and cross-correlations. The advantage of this test over 
reconstructed Cf'^ is that we are able to see tomographic 
view of consistency relation. 



V. DISCUSSION 

In this paper, we advocate the use of consistency re- 
lations to test gravity on cosmological scales. We did so 
by using a combination of obscrvables. This approach 
remains model independent and does not rely on any 
specific parametrization. We focused in particular on 
the joint use of galaxy surveys and weak-lensing observ- 
ables. We showed how using the former to predict the 
latter, we can build a strong self-consistency test for GR. 
We also showed how large-scale redshift surveys can be 
extremely valuable in such an endeavor. The test we 
proposed seems particularly appealing since any weak- 
lensing survey is also by nature a galaxy clustering survey 
and we illustrated how even a photo-z survey can be used 
to build a strong self-consistency test. We thus expect 
this test to be applied in the near future to the CFHTLS 
survey ^ and others hke DES the proposed JDEM 
[sHi and Euclid [6^ space missions. 

We address in this work the key observational and as- 
trophysical systematics like the galaxy bias and redshift 
uncertainties but many more survey specific ones would 
have to be carefully studied before any conclusion can 
be drawn from any particular data-sets, e.g. n{z) errors, 
PSF correction and variation over large scales. It is quite 
likely that this kind of test will be ultimately limited by 
the level of control of systematics. In particular, we pur- 
posefully worked on linear scales. Non-linear effects in 
real and redshift space will limit the useful scales for such 
a program. If we assume that we can handle non-linear 
effects up to k=0.3 h/Mpc, we can see in the left panel of 
Fig. m that our program is still tractable. Besides, weak- 
gravitational lensing has now been measured in the linear 
regime where this test could be applied [ss'] . Another un- 
explored sides of our work lie in the possible degeneracies 
with underlying cosmological model parameters but we 
leave this study for future work. 

More positively, although, our survey parametrization 
is quite generic, these surveys could certainly be cus- 
tomized to enhance the discriminatory power of this test 
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and others. We encourage future missions to include the 
feasibihty of such tests as an extra criteria in the opti- 
mization of their design. 

Assuming that aU systematics are weh under control, 
we can wonder what new physics can be probed with 
this test. As an illustration, we plotted in Fig. [6]the pre- 
dictions for a f{R) theory where a non-linear function 
/ of the Ricci scalar R is added to the Einstein- Hilbert 
action [s^. We considered typical value of the dimen- 
sionless Compton wavelength parameter Bq — 1 (see Eq. 
17 in [5^). In /(i?) gravity models, the Poisson equation 
is modified so that Gn ^ Gm/{^ + df/dR) and a non- 
zero anisotropy stress is introduced to break the simplest 
GR model anisotropy condition <i> = — 5*. As visible in 
Fig- El the predicted Cf''' for a f{R) gravity model using 
GR assumptions would differ from the observed one at 
a level well exceeding the various potential biases. This 
discrepancy originates from the GR assumption in Eq. l34l 
since any non-zero anisotropic stress invalidates the re- 
lation 2$ = $ — '5. However, the reduced mass scale 
affecting f{R) theories modify the Poisson equation. It 
can thus be probed directly using the cross-correlation 
technique detailled in Sec. lIVBl Note that since the dif- 
ference between this f{R) model and GR increases when 
the redshift decreases, so does the discrepancy between 
the observed and predicted Cf'^ for a f{R) theory. From 
the Fig. [7] it can be seen that constraints on Bq of order 
1 or less are in principle feasible with such a test. Other 
aspects of f{R) theories are discussed in [stI IssI. [59l|. 



More generally, it can be argued that any theory with 
non-minimal interaction in the dark sector or dark energy 
clumping would fail this test since it either modifies the 
Poisson equation or introduce an anisotropic stress (3^ . 
iGOi] . Reversely, it is important to note that some types of 
modified gravity models would pass the test we detailed 
here. For example, the potentials following from a DGP 
model [6lj are modified in such a way that the light path 
tracing potential differences is similar to the standard 
ACDM predictions. To detect such models, we need to 
implement the program proposed in [4^ . 

The tests presented above are only two of examples of 
the possible one highlighted in Sec. [Til and we plan to 
extend our studies to include other consistency relations 
as well as other observables like the integrated Sachs- 
Wolfe effect 62| and the lensing of the cosmic microwave 
background. We believe however that the observational 
program sketched in this paper already offers a stringent 
observational test of GR and we thus plan to put it into 
action using current data. 
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